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■ We consider a small and fixed number of fermions in an isolated one- 
\Q • dimensional trap (microcanonical ensemble). The ground state of the sys- 
tem is defined at T = 0, with the lowest single-particle levels occupied. The 
number of particles in this ground state fluctuates as a function of excitation 
energy. By breaking up the energy spectrum into particle and hole sectors, 

■ and mapping the problem onto the classic number partitioning theory, we for- 
mulate a new method to calculate the exact particle number fluctuation more 
efficiently than the direct combinatorics method. The exact ground state 

fS| ! number fluctuation for particles interacting via an inverse-square pair-wise 

interaction is also calculated. 
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I. INTRODUCTION 



Recent experimental observation of quantum degeneracy in a dilute gas of trapped 
fermionic atoms at low temperatures by DeMarco and Jin |I[ inspired our previous work 
on fermions ||. We considered an ideal gas of spinless fermions in an isolated harmonic 
trap in one- and two-dimensions. In this work we restrict our study to a system of fermi 
gas in one- dimension. In addition to reporting some new results we obtained for a system 
of ideal fermi gas, we now introduce inter-particle interaction. Since our system is isolated 
(microcanonical ensemble), the results given in this paper for both ideal and interacting 
fermi gases are exact. 

Consider first a system containing N non-interacting particles in a one-dimensional har- 
monic confinement, either bosons or (spinless) fermions. The canonical partition function 
of this system is given by 

N i oo 
Z N (x) = X E «® T] = X E N (0) £ n(nj N)x n^ (1) 

j=l I 1 _ XJ ) n=0 

where x = e _/3 , (3 = 1/ksT, -Ejv(O) is the ground state energy of the system. En(0) = N/2 
for bosons and E^(0) = N 2 /2 for fermions. For simplicity we have set tkv = 1, where uj is 
the oscillator frequency. Note that the expansion coefficient Q(n, N) is the same for bosons 
and fermions, and so is the canonical partition function, apart from the prefactor involving 
the ground state energy E N (0). The ground state of the system is defined at T = 0. For 
fermions the ground state consists of a set of N occupied energy levels, while for bosons it is 
the single lowest energy level. Given an excitation energy n, there are many ways in which 
the particles can share this energy and excite to the higher states, so that the particle number 
in the ground state fluctuates. This ground state number fluctuation for ideal particles may 
be calculated exactly using combinatorics. This was done by us in Ref. for ideal fermions. 
The desired quantity is the microcanonical multiplicity u(n, N ex , N), defined as the number 
of ways to distribute exactly N ex out of N particles from the ground state to the excited 
states, given n quanta of energy in a system of N particles. Once lu(ti, N ex , N) is found, the 
corresponding probability P(n, N ex , N) and the exact ground state number fluctuation can 
be calculated straightforwardly: 

P { n,N^N)= (2) 
E Nex =iu{n } N eX} N) 

N 

(N ex ) = J2 P(n,N ex ,N)N ex , (3) 

N em =l 
N 

(N 2 ex )= £ P(n,N ex ,N)Ni, (4) 

N ex =l 

m f = {nd - (au 2 

= (N 2 ) - (iV ) 2 . (5) 

Note that by definition P(n, N ex , N) is normalized to unity. Eqs. (@)-(|5p apply to any system 
of particles, either bosons or fermions. Note, however, that the microcanonical multiplicity 
u(n, N ex , N) is different for different types of particles. For fermions, apart from the fact that 
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the ground state consists of more than one energy levels, the distribution of N ex particles 
must comply with the Pauli Principle. For clarity, therefore, we shall attach a bracketed 
superscript (n, N ex , N) for bosons, and uj^ F ^{n,N ex ,N) for fermions in what follows. 

Although the combinatorial method may be used to determine the microcanonical mul- 
tiplicities u B (n, N ex , N) and u) F {n, N ex , N) for both ideal bose and fermi gases exactly, it is 
very time-consuming computationally. In the case of an ideal boson gas in a one-dimensional 
harmonic trap, the problem is greatly simplified due to its connection to number partitioning 
theory. This is because the energy spectrum of a one- dimensional harmonic trap is equally 
spaced, and the N ex bosons are excited from a single lowest energy level. Given n excitation 
quanta imparted to the system, the number of ways to excite exactly N ex bosons to the 
excited states, ui B {n, N ex , N), is equivalent to the number of ways of partitioning an integer 
n into exactly N ex parts. It is related to the expansion coefficient Q(n, N) of the canonical 
partition function Z^{x) via 0: 

u B {n, N ex , N) = Q{n - N ex , N ex ), n > N ex 

= 0, otherwise. (6) 

Note that Eq. (||) implies that u{n, N ex , N) = u{n, N ex , N'), N' > N ex , i.e. independent of 
system size. In the language of number theory, the canonical multiplicity Q(n, N) is the 
number of ways of partitioning n into 1, 2,3, ...,N parts. Clearly, 

TV 

n(n,N) = £ u B (n,N ex ,N). (7) 

N ex =l 

Since the canonical multiplicity Q(n,N) is the same for bosons and fermions (Eq. ([[])), the 
following must be true: 

N 

n(n,N) = ]T co F (n,N ex ,N). (8) 

N ex =l 

So even though the bosonic uj B {n, N ex , N) is very different from the fermionic u) F {n, N ex , N), 
their sum Q(n, N) is the same. For bosons, expression @ implies that knowing the canonical 
multiplicity Q(n, N) is sufficient to study the system at a microcanonical level. It therefore 
allows an easy way to calculate the exact number fluctuation, instead of the more cumber- 
some combinatorics formulae. This is not the case for fermions, however, since Eq. ([]) does 
not apply. This is because the (quantum mechanical) ground state of fermions consists of 
occupied energy levels, unlike bosons whose ground state consists of only a single lowest 
energy level. To see why this multi-level ground state of fermionic system prevents the 
application of number theory, consider a system containing N = 2 particles. For simplicity 
take n = 2 and consider exciting one particle, i.e., N ex = 1. In the language of number 
theory, this means partitioning the integer n = 2 into one part. There is only one way of 
doing this. For bosons there is also only one way of taking 1 out of = 2 particles from 
the ground state, and put it in the second level above, i.e., u> B (2, 1,2) = 1. For fermions, 
however, either the particle which lies in the fermi level Ep may get excited to the second 
level above, or the particle which lies just below Ep gets excited to the first level above. 
Thus uj f (2, 1, 2) = 2, which is not equal to one as in the case of bosons. Note that n = 2 may 
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also be partitioned into two parts, and Q(2, 2) = 2 since 2 = 2, 1 + 1. For bosons this means 
exciting both particles, i.e. u B (2, 2, 2) = 1, whereas for fermions u F (2, 2, 2) = since it takes 
at least 4 quanta to excite 2 fermions. In either case the microcanonical multiplicities sum 
up to 2 (Eqs. (^)and (§)). So the multi- level ground state of fermionic system compensates 
for the restriction in the number of ways the particles may be distributed imposed by the 
Pauli Principle. Fortunately, the mapping of the system of N fermions in a one-dimensional 
harmonic trap to number partitioning theory is not impossible. In section [TJ, we show that 
by looking at the excited energy levels and the multi-level ground state separately, we are 
able to connect each of these two spaces to a similar bosonic problem, which in turn pertains 
to number partitioning. We successfully derive a formula similar Eq. (|6|) which applies for 
fermions. Apart from mathematical interest, the formula allows us to calculate the exact 
ground state number fluctuation more efficiently than the combinatorics method. 

In the conventional treatment of the ideal bose gas, the number fluctuation is calculated 
using the grand canonical ensemble (GCE). The result is given by M 



E K) (i + K» 

. k 



1/2 

(9) 



where (rife) is the occupancy of level k. As the temperature approaches zero, the particles 
condense to the ground state. Relation @ thus yields ((5N )) rs (Nq) ~ (N). This result 
is unphysical since as T — > 0, there is no energy to excite particles and hence the number 
fluctuation should vanish. This is a failure of the GCE in describing the fluctuation of an 
ideal bose gas at low temperatures. Experimentally, the system has a fixed total number 
of particles N confined in space by a trapping potential, so that it is more accurately 
described by a canonical ensemble (CE), or better still, a microcanonical ensemble (MCE). 
This has been shown by several authors; some used property (||) for the MCE, and the 
ground state number fluctuation has been found to vanish as T — > ||. Within the GCE 
framework, however, the divergence can be removed by introducing inter-particle interactions 
||. Consider an ideal gas in a box, and recall that in GCE the particle number fluctuation 
is related to the thermal compressibility via 

1 (dV\ V (5N\ 2 , . 

For T —>■ 0, the pressure is independent of volume, so that \ T —>■ oo and hence 5N — > oo. 
With a repulsive interaction, however weak, there is a pressure due to interaction to ensure 
the fluctuation to be finite. In the case of fermions, there exists a Pauli pressure even if the 
gas is ideal so that the ground state fluctuation finite [ffl: 



(SN t 



o 



nl/2 



(11) 



the sum k runs over the ground state levels occupied at T = 0. As T —>■ 0, (rife) — > 1 
for k below fermi level, so that the fluctuation goes to zero as expected. In this paper we 
shall use the inverse square two-body interaction in one-dimension (the Calogero-Sutherland 
model, abbreviated CSM) . This interacting model may be mapped onto the Haldane-Wu 
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generalized exclusion statistics (also known as fractional exclusion statistics, FES) whose 
quasiparticles are non-interacting, characterized by a parameter g ||8HTC|1. g = 0,1 corre- 
spond to bosons, fermions respectively. Other values of g represent various degrees of 'Pauli 
blocking'. The model is solvable, and the energy spectrum is exactly known. In |TT| we 
calculated the fluctuation for bosons and fermions using this interaction in both GCE and 
CE. For an interacting bose gas we showed that the divergence of the number fluctuation 
at T = in the GCE is removed. The ground state number fluctuations of an interacting 
gas using CE in two- and three-dimensions have been previously calculated by a number of 



authors for different types of interaction [TT,TI|. Our choice of the inverse-square two-body 
interaction model in one-dimension has a two-fold advantage. First, it serves as an interact- 
ing model which allows us to demonstrate that the divergence in the GCE is removed when 
an interaction is introduced. Second, owing to its mapping to non-interacting FES particles, 
the model enables us to show that the divergence is removed even for an ideal gas provided 
that the Pauli blocking is not zero, no matter how weak. To the best of our knowledge, the 
exact fluctuation of interacting particles in the MCE has never been calculated before. In 
section |T| of this paper, we report some results of the exact fluctuation of trapped particles 
interacting via CSM, or equivalently, of ideal FES particles. 



II. IDEAL FERMIONS 

Consider N spinless fermions in a one-dimensional harmonic confinement. The ground 
state is defined at T = where the N particles occupy the lowest energy levels up to and 
including the fermi level Ep. Given n quanta of energy, consider breaking up n into two 
parts: 

n = n h + n p , (12) 

where is the number of quanta it takes to bring N ex particles to the fermi level Ep (which 
is equivalent to the distribution of N ex holes to the states below and including Ep), and n p 
is the number of quanta it takes to distribute these N ex particles in the excited states (above 
Ep). This effectively divides the fermionic energy levels into two sectors: the particle space 
above E F and the hole space below and including E F . Let u^(n h , N ex , N) be the number 
of ways to distribute N ex holes in the hole space, and u^(n p , N ex , N) be the number of 
ways to distribute N ex particles in the particle space, both according to Pauli Principle, then 
given n quanta and N ex particles, 

u F (n,N ex ,N) = J2 oo { h F} (n h ,N ex ,N)uj { p F} (n p ,N ex ,N), (13) 

where the set {nh,n p } satisfies Eq. (|i~2f) for a given n. The problem now pertains to finding 
(jjf^foh, N ex , N) and uo p {^n p ,N ex ,N). At first glance this seems to be more complicated 
than finding a single quantity u F (n, N ex , N). However, recall that u F (n, N ex , N) is the 
number of ways to distribute N ex particles above Ep with respect to a set of N ground 
state energy levels. By breaking up the fermionic energy levels into two parts we are now 
distributing N ex particles above Ep and N ex holes below Ep, both with respect to a single 
energy level. As we shall see shortly, this allows us to connect to number partitioning theory 
as in the case of bosons. Let us now look at these two spaces separately. 
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A. Particle space 



First consider the particle space. This space is unbounded starting from the fermi level 
Ep. We now consider a situation in which the N ex particles have already been taken out 
of the hole space, and are residing in the fermi level Ep awaiting to be distributed in the 
particle space. Note that n p quanta are distributed among N ex fermions, with respect to only 
one energy level Ep. The problem is now similar to the bosonic one, except the distribution 
of particles must comply with the Pauli Principle. In terms of number partitioning theory, 
N ex , N) is the number of ways of partitioning an integer n p into N ex different parts, 



with N ex < N. As in the case of bosons, it is related to the canonical multiplicity via ||: 

uj { /\n p , N ex , N) = n(n p - A p , N ex ), 

= n(n p K — L ,N ex ). (14) 

Remarkably, this is the same form as Eq. (|6]) for bosons with the shifted energy N ex replaced 
by A p . This shifted energy, A p = N ex (N ex + l)/2, is in fact the minimum energy it takes to 
excite N ex particles from Ep, 

n™ n (N ex ) = (iVe * + 2 1)iVex . (15) 

Note that the partition function of this space is no longer given by that of fermions. The 
canonical multiplicity in Eq. fll4"| ) should therefore be thought of as the bosonic multiplic- 
ity. In other words, the problem is now mapped onto a similar bosonic problem, with the 
restriction that the parts of an integer being partitioned are different. This notion is most 
helpful when we discuss the hole space. 



B. Hole space 

We now consider taking N ex particles out of the multi-level ground state and put them 
in the fermi level Ep (or equivalently, creating N ex holes in the ground state). For a given 
n h quanta, we wish to find uo^(n h , N ex , N), the number of ways of doing this. Unlike the 
particle space, the dimension of the hole space is bounded, set by the number of particle N 
of the system. Note that for a given number N ex of particles, the Hilbert space of available 
states for N ex holes is dependent on the value of N ex itself and is given by: 

N H = N-N ex . (16) 

We need to find the partition function of this space for each N ex , and derive a formula 
similar to Eq. flT4]) for ujjf^ (n,h, N ex , N). This may be done by considering a new system 
containing N ex bosons, whose energy space is bounded and is given by Nh + 1 including the 
ground state. The goal is to determine Z^ ex (x), and expand this in terms of the coefficient 
Qh{rih, N ex ). The A ex -hole partition function of this hypothetical system can be found using 

i N ex 

zU*) = j^Y,z?U*)zL-M ( 17 ) 



6 



where Z^(x) is the one-particle partition function of the system containing N ex bosons and 
is given by: 

N„ 

Zf(s)=5>\ (18) 

i=0 

Note that the one-particle partition function needs to be determined for a given N ex . Using 
this, Z§ ex (x) may then be found from ([[7]). Once Z^ ex (x) is found, we may expand it in 
powers of x: 



7 u 



[X 



J2tt h (i,N ex )x\ (19) 



We are now ready to determine a formula for u\ (nh, N ex , N). Because the hole space 
includes Ep, the minimum energy to create a hole (or dig a particle) to put in the fermi 
level is zero, since there already is a particle there; for two holes the minimum energy is one, 
for three holes it is three... etc. In general, 

nT n {N ex ) = (jkzWgL . (20) 

Similar to Eq. ( |14|) with the energy shift Ah given by (0), the number of ways of creating 
N ex holes in the hole space oj^ijih, N ex , N) is given by: 

co h (n h , N ex , N) = n h (n h - N ^ N ™~ 1 ) ^ ^ {2l) 

Using Eqs. (|14]) and (0), the number of ways of distributing N ex fermions to the excited 
states, Eq. (|T3|), now reads: 



co F (n,N ex ,N)= £ 0(n p ~ Nex{N ™ + 1} , N e «) ~ ^ ^ , N ex ). (22) 

It is obvious that if there is no hole space, u>h(n h , N ex , N) = Qh( n h — A h , N ex ) = 1, the 
sum over the set {ri/j, n p } vanishes since n = n p , the energy shift A p = nT ,%n (N ex ) = N ex , 
and Eq. (|13|) reduces to Eq. (|) for bosons. 

Using Eqs. (^) and (|2])-(|5]), we calculated the ground state number fluctuation of 
fermions for N = 100. The result is shown in Fig. p]. For comparison we also show the 
corresponding result in the CE. 



III. IDEAL FES PARTICLES 



We now calculate the exact ground state fluctuation of particles interacting via CSM, or 
equivalently, non-interacting particles obeying FES [fFHIOfl. The Hamiltonian of the CSM is 
given by 0: 



TV 



i=l 



h 2 d 2 l 2 , 

2^ + 2 m " X * 



h 2 N 

+ - E 



m yr . [Xi 

1<J = 1 v ' 



(23) 
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with the dimensionless coupling parameter g > 0. The particles are confined in a harmonic 
well and the thermodynamic limit is obtained by taking u — > as N — > oo, with uoN = 
constant. In the thermodynamic limit, the properties of the system are translationally 
invariant, and would be the same if the particles were on a line, or a circle, instead of a 
harmonic confinement. To make the problem well-defined quantum mechanically, we have 
to demand that the wave functions go to zero as \x{ — Xj\ 9 whenever two particles i and j 
approach each other. Since the particles cannot cross each other, we may choose the wave 
function to be either symmetric (bosonic) or antisymmetric (fermionic). For g = and 1, 
the model describes free bosons and free fermions respectively. 

In [[[I] we have calculated the fluctuation for FES particles in the GCE and CE. To do 
this in the MCE, we need to find cu^(n, N ex , N) for a given g. The energy spectrum of the 
CSM Hamiltonian (p3|) is exactly known, and may be expressed in terms of the quasiparticle 



energy spectrum e g of non-interacting FES particles, which is given by §J: 

e g = e k -(l-g)N k , (24) 

with Hu = 1. Nk is defined as the number of particles below energy level k, = (k — 
|),A; = l,2,...is the harmonic oscillator energy spectrum. The canonical N-particle partition 
function of FES particles reads: 

z N = x 9} (o) n 7 -i 

j=i (1 

oo 

= x E » W £ n(n, N)x n . (25) 

n=0 

Note that the canonical multiplicity Q(n, N) is still the same for FES particles as for bosons 
and fermions. The only effect of the interaction strength g is to alter the overall ground 
state energy E^(0) = gN(N — l)/2 + N/2 which reduces to that of fermions, bosons for 
(7 = 1,0 respectively. 

Using Eq. (p4j) , one may determine the spectra of FES particles for any value of g. In the 
Fig. (Q) we show an example of how to find the spectrum of semions (g — 1/2) for some values 
of excitation quanta n and number of particles N = 5. Starting from the fermionic spectrum 
on the left, for each value of n, the spectrum of semions is drawn and the resulting particles 
which are in the excited states are determined. In the case of n = 3, for instance, the only 
possibility for fermions is to excite one particle and there are 3 ways of doing this as shown. 
Therefore, cu^(3, 1,5) = 3, and c</ F }(3, 2,5) = to^(3, 3, 5) = since 3 quanta is too few 
to excite 2 or more particles. For semions, however, u^ 1 / 2 '^, 1,5) = 1, 10^^^(3,2, 5) = 2, 
and a;' 1 ' 2 } (3, 3, 5) = 0. Note that fi(3, 5) = 3 in both cases. Following the same pro- 
cedure, we found the microcanonical multiplicities u^(n,N ex ,N) for g = 3/4,1/2,1/4, 
n = 1...16 quanta, and iV = 2. ..5 particles. The method gets more cumbersome, however, 
for larger values of n and N. Ideally, one wishes to be able to determine 
from co^ F ^(n, N ex , N) or to^ B ^(n, N ex , N), or from £l{n,N ex ) similar to Eqs. @ and (p2|) . 
This general formula for u^(n, N ex , N) for any g, if it exists, is yet to be found. Here, we 
report the finding of the microcanonical multiplicities for only two values of g: g = 

(close to fermions), and g = jt^t (close to bosons). For g = j^zj, a/' ?rrT ^(n, N ex , N) is found 
to be given by: 



S 



For g 



oo { ^ } (n, N ex , N) = wW(n + N, N ex , N) - u^ F \n + N, N ex , N - 1), (26) 
1 . 



N-V 



u { ^\n, N ex , N) = u {B} (n, N ex , N), N ex ^N-1,N 

uj { ^ } (n, N - 1, N) = u {B} (n, N - 1, N) + u {B} (n - 1, N - 1, N), 
uj { ^ } (n,N,N) =uj {B} (n- N,N,N). (27) 

These multiplicities must, of course, satisfy 

N 

tt(n,N) = ]T u^(n,N ex ,N). (28) 

N ex =l 

For < 5 and n < 16, the values found using Eq. (p6|)-(p7|) were verified with those found 
using the direct combinatorial method as described above. 

Although a general formula for the microcanonical multiplicities for any g has not been 
found, an important point concerning them is observed. For a given number of particles N, 
consider a set of discrete values of g, 

N-2 N-3 1 , s 

s = 1 —v—v-—v°- (29) 

Since the levels of the FES particles are shifted by an amount which depends on the value of g 
and the number of particles below (see Eq. (j2"3D), for some values of g a particle might lie very 



close to the last level of the ground state, Ep . However, the particles are considered excited 
if and only if they lie above Ep 9 ^, no matter how close. This results in the multiplicities 
cu^(n, N ex , N), for > g > , i — 1, N — 1 to be the same. Note the equal sign 

in 9 > N at-1 ■ F° r instance, 



uo {9 \n,N ex ,N) = oj^(n,N ex ,N), for > g > 0. 

Note that for these values of g given by Eq. (||), the 'Fermi level' E { F 9} = N- (1 -g)(N- 1) 
is integral. So for g = (N — k)/(N — 1), where k is an integer, the 'Fermi level' lines up 
with the (k — l) th level below Ep of fermions. In box 1, n — of Fig. for instance, where 
N = 5 and g = 2/4 (k = 3), the 'Fermi level' lines up with the second level below Ep. This 
explains why the microcanonical multiplicities for some range of g such that the 'Fermi level' 
lies between the k th and (k + l) th levels of fermions are the same. This is due the discrete 
nature of the energy levels. 



IV. DISCUSSION 

In our previous work (Ref. [[J) where the direct combinatorial method was used, the 
fermionic calculation of the ground state fluctuation was restricted to a low number of 
particles and quanta n using a normal office computer (Pentium III, 500 cpu). For a 
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relatively small number of particles (e.g., N = 10), at higher excitation n the combinatorics 
method is more time-consuming due to the rapid increase in the number of possibilities 
with n. The method described in section f[I| translates the problem in combinatorics into 
the problem in calculating the partition functions of the hole space, the latter being simpler 
computationally. Although this method is still time-consuming and the calculation for larger 
N (N > 100) is still not possible using our office computer, it is more effective for higher 
number of quanta and relatively small number of particles. For demonstration we display 
the ground state fluctuation of fermions as a function of energy quanta n in Fig. [I] for 
N = 30. We also show the corresponding curve in CE, which is the same as GCE except 
at very low temperature (see e.g., Ref. 0) for comparison. As expected, both go to zero 
as T — > 0, with the MCE fluctuation less than that given by CE for all n. Note that the 
two fluctuations are very different, even for very high excitations. At n = 6000 which is 
200 x Ep, the CE curve still differs from the MCE one by about 14%. It was shown in 
Ref. 0] for iV = 15 that the ground state occupancy (N Q ) = J2k F ( n k) for the two ensembles 
are very similar. Clearly, the number fluctuation is more sensitive to the ensemble used. 
Therefore, for a relatively small particle number, while it may be adequate to use the CE (or 
GCE) to describe a thermodynamic quantity such as the ground state occupancy, it should 
be used with caution when calculating the number fluctuation and related quantities. 

Using Eqs. (|26|), QSTD and (@)-(0), we calculated the exact ground state number fluctuation 
for interacting particles. In Fig. |3]a we display the fluctuation for N = 5, g = 3/4,1/4. 
For comparison we also show the fluctuation of free fermions and bosons. The values for 
g = 3/4, 1/4 obtained using the direct combinatorial method (see Fig. ||D for n < 16 are also 
shown. Note that the curves agree with these values exactly. Note also that the number 
fluctuation for free fermions and free bosons cross at a certain energy, with the fermionic one 
starting from smaller at small quanta, to larger at high quanta. This is because the number 
of possibilities of creating holes within the fermi sea and distributing particles above, which 
starts from low at small energy, increases more rapidly than for bosons whose ground state 
consists of only one level. Similar behaviours are observed for FES particles, whose g values 
(g — 3/4, 1/4) represent partial Pauli blocking which are both less than that of fermions. 

Eqs. fl26|) and ( p?|) in principle may be applied for any N. However, since they involve the 
addition of two quantities which may be very large at large quanta and number of particles, 
there is a difficulty in obtaining their values accurately. Therefore, without lost of accuracy, 
we restricted the calculations to iV = 10. Fig. shows the ground state number fluctuation 
for N = 10, g = 1, 8/9, 1/9, 0. Note that the curves for g = 8/9, 1/9 are closer to those of 
fermions and bosons than for g = 3/4, 1/4, N = 5, in Fig. |3|a. As N gets larger we expect 
the fluctuation curves for g = ^=1 and jr-r to come very close to those of free femions and 
free bosons. The formulae (|26|) and (|27|) are therefore useful only for systems with a small 
number of particles. Note also that in both graphs, the results for g = are closer to those 
of fermions than for g = -^-j- to bosons. This is can be understood from comparing Eqs. ( |26|) 
and (pTp. Eq. (^) involves the difference of the fermionic microcanonical multiplicities of 
two system sizes, whose values may be very similar especially at low quanta. Eq. (p7|), 
however, involves the addition of two bosonic microcanonical multiplicities. This brings the 
microcanonical multiplicities of g = -^-j- further from those of bosons than g = horn 
those of fermions. 

Finally, In Figs. f|a,b we compare these exact fluctuation curves with the corresponding 
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CE curves for g = ^zf > ]vZT> ^ = 5, 10. From fll |, the expression for the canonical ground 
state fluctuation is given by: 

((SN ) 2 ) {9} = g((5N )Y + (1 - g)((5N ) 2 ) B (30) 

Clearly, the results between the two ensembles are very different especially for small n. 
In both graphs the microcanonical fluctuations are less than the canonical ones for all n. 
Note that the CE formula (|30| ) allows one to calculate the fluctuation for a general g. The 
determination of a general formula for the microcanonical multiplicity, to {9} (n,N ex ,N), and 
hence the ground state number fluctuation for any g, if it exists, remains a challenge. 
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FIGURES 



FIG. 1. Ground state number fluctuation of fermions as a function of excitation energy n (in 
unit of huS) for N = 30. The result in the CE is also shown for comparison. 

FIG. 2. Spectra of semions (g = 1/2) derived from those of fermions using Eq. (|24|). The 
(quantum mechanical) ground state is defined at n = 0. The spectra for different energies are 
separated by boxes. Within a box there may be more than one configuration, each of which 
contains the spectrum of ideal fermions on the left and semions on the right. The number of 
excited particles, N ex , are indicated in each case. 

FIG. 3. (a) Ground state number fluctuation curves as functions of excitation quanta n for 
N = 5, g = 1,3/4,1/4,0. For clarity the low energy part is shown in the inset. The data 
represented by the symbols are obtained using the direct combinatorial method as shown in Fig. ^; 
O for g = 3/4, and □ for g = 1/4. (b) Same as in (a), for N = 10 and g = 1,8/9,1/9,0. 

FIG. 4. (a) Comparison between the exact fluctuations as discussed in this paper with those in 



the CE (Eq. (30)) for N = 5, g = 3/4, 1/4. (b) Same as in (a), N = 10, g = 8/9, 1/9. 
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